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Two  objectives  dominate  consideration  of  control  moment  gyroscopes  for  spacecraft  maneuvers:  high  torque 
(equivalently  momentum )  and  singularity-free  operations.  This  paper  adds  to  the  significant  body  of  research  toward 
these  two  goals  using  a  minimal  three-control-moment-gyroscope  array  to  provide  significant  singularity-free 
momentum  performance  increase  spherically  (in  all  directions)  by  modification  of  control-moment-gyroscope  skew 
angles,  compared  with  the  ubiquitous  pyramid  geometry  skewed  at  54.73  deg.  Spherical  1 H  (one  control  moment 
gyroscope’s  worth  momentum)  singularity-free  momentum  is  established  with  bidirectional  1 H  and  2 H  in  the  third 
direction  in  a  baseline  configuration.  Next,  momentum  space  reshaping  is  shown  via  mixed  skew  angles  permitting 
orientation  of  maximum  singularity-free  angular  momentum  into  the  desired  direction  of  maneuver  (yaw  in  this 
study).  Finally,  a  decoupled  gimbal  angle  calculation  technique  is  shown  to  avoid  loss  of  attitude  control  associated 
with  singular  matrix  inversion.  This  technique  permits  3 H  (maximal )  y a w  maneuvers  without  loss  of  attitude  control 
despite  passing  through  singularity.  These  claims  are  demonstrated  analytically,  then  heuristieally,  and  finally 
validated  experimentally. 
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Nomenclature 

gradient  matrix  of  gimbal  angles  and  skew  angle(s) 
with  respect  to  gimbal  rotation  angle 
cosine  function 
determinant 

total  momentum  vector  {  hx  hy  hz  }T 
magnitude  of  total  angular  momentum  of  the  array 
of  control  moment  gyroscopes 
total  momentum  magnitude  of  control-moment- 
gyroscope  array  when  encountering  a  singular 
condition 

momentum  in  the  x,  y,  and  z  directions  (normalized 
by  one  control  moment  gyroscope’s  worth  of 
momentum) 
sine  function 

skew  angle  for  ith  control  moment  gyroscope 
Partial  derivative  of  momentum  in  x-direction  taken 
with  respect  to  ith  gimbal  rotation  angle 
gimbal  rotation  angle  for  ith  control  moment 
gyroscope 

gimbal  rotation  rate  for  ith  control  moment 
gyroscope 

normalized  angular  momentum  output  by  one,  two, 
and  three  control  moment  gyroscopes 
normalized  angular  momentum  output  by  four 
control  moment  gyroscopes  (also  called  saturation) 


I.  Introduction 

RAPID  spacecraft  reorientation  often  drives  design  engineers  to 
consider  control  moment  gyroscopes  (CMGs).  CMGs  are 
momentum  exchange  devices  that  exhibit  extreme  torque  mag¬ 
nification  (i.e.,  for  a  small  amount  of  torque  input  to  the  CMG  gimbal 
motors,  a  large  resultant  output  torque  is  achieved)  but  inherently 
possess  singular  directions,  where  no  torque  can  be  generated. 
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Despite  singularity  issues,  CMG  research  began  in  the  1960s  for 
large  satellites  like  Skylab  (which  used  three  double-gimbaled 
CMGs,  or  DGCMGs)  [1],  Computers  of  the  time  could  not  perform 
matrix  inversion  in  real  time.  Simple  systems  that  did  not  require 
matrix  inversion  were  an  obvious  choice.  Otherwise,  algo¬ 
rithmically  simple  approximations  must  have  been  available  for 
the  system  chosen.  CMG  steering  and  singularity  avoidance  was 
researched  a  lot  in  the  1970s,  1980s,  and  1990s  [2-11], 

Singularity  avoidance  was  typically  done  using  a  gradient  method 
and  DGCMGs  [5,6rf2].  These  gradient  methods  are  not  effective  for 
single-gimbaled  CMGs  (SGCMGs)  like  they  are  for  DGCMGs. 
Magulies  and  Aubrun  were  the  first  to  formulate  a  theory  of 
singularity  and  control  [7]  including  the  geometric  theory  of  singular 
surfaces,  the  generalized  solution  of  the  output  equation,  null  motion 
(using  greater  than  three  SGCMGs),  and  the  possibility  of  singularity 
avoidance  for  general  SGCMG  systems.  Also,  in  1978,  Russian 
researcher  Tokar  published  singularity  surface  shape  description, 
size  of  work  space,  and  considerations  of  gimbal  limits  [8], 

Kurokawa  et  al.  identified  that  a  system  such  as  a  pyramid  type 
SGCMG  system  will  contain  an  impassable  singular  surface  and 
concluded  systems  with  no  less  than  six  units  provide  adequate  work 
space  free  of  impassable  singular  surfaces  [17] .  The  Mir  space  station 
was  designed  for  six-SGCMG  operations.  The  research  contained  in 
this  paper  evaluates  singularity-free  operations  using  a  mere  three- 
CMG  array  by  reducing  singularities  and  then  penetrating  those 
remaining  singularities  without  loss  of  attitude  control. 

Continued  research  aimed  at  improving  results  with  less  than  six 
CMGs  emphasized  a  four-CMG  pyramid  [  14].  Much  research 
resulted  in  gradient  methods  that  regarded  passability  as  a  local 
problem  that  proved  problematic  [l_3,j_5,lA].  Global  optimization 
was  also  attempted  but  proved  problematic  in  computer  simulations 
[17-19].  Difficulties  in  global  steering  were  also  revealed  in  [20], 
Reference  [21]  compared  six  different  independently  developed 
steering  laws  for  pyramid-type  SGCMG  systems.  The  study 
concluded  that  exact  inverse  calculation  was  necessary  (the  exact 
inverse  calculation  is  used  in  this  paper).  Other  researchers  addressed 
the  inverted  matrix  itself,  adding  components  that  made  the  matrix 
robust  to  inversion  singularity  [22-24]  as  extensions  of  the  approach 
to  minimize  the  error  in  generalized  inverse  Jacobian  calculation 
[25].  These  approaches  introduced  tracking  errors  where  necessary 
to  avoid  singularities  by  following  a  different  momentum  path 
(generating  other-than-desired  torque).  Reference  [26]  sought  to  use 
a  hybrid  steering  logic  to  maintain  attitude  tracking  precision  while 
avoiding  hyperbolic  internal  singularities  or  escaping  elliptic 
singularities  with  a  four-CMG  array  skewed  at  54.73  deg,  the 
ubiquitous  baseline  described  in  [141  and  depicted  in  Fig.  1 . 
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Fig.  1  Singularity  surfaces  associated  with  the  ubiquitous  skew  angle 
of  54.73  deg. 


succinctly  stated:  1)  optimize  spherical  momentum,  and  then 
2)  minimize  impact  of  singularities.  The  approach  adopted  here  will 
reverse  the  traditional  approach  as  follows:  1 )  minimize  singularities, 
and  then  2)  maximize  spherical  momentum.  Surprisingly,  the  result 
turns  out  quite  differently.  This  paper  begins  by  comparing  perfor¬ 
mance  of  a  three  of  four  (3/4)-skewed  pyramid  using  the  baseline 
skew  angle  /l  =  54.73  deg  (derived  from  the  four-CMG  pyramid)  to 
the  performance  of  a  3/4  CMG  pyramid  with  the  skew  angle 
specifically  optimized  for  maximum  singularity-free  angular 
momentum.  It  will  be  shown  that  a  3/4  CMG  pyramid  skewed  at 
/I  =  54.73  deg  can  achieve  a  maximum  singularity-free  momentum 
of  15%  of  one  CMG’s  maximal  momentum,  while  the  proposed 
geometry  can  achieve  100%  of  one  CMG’s  maximal  momentum  in 
all  directions  with  superior  performance  in  a  preferred  direction. 


Momentum  path  planning  is  another  approach  used  to  attempt  to 
avoid  singularities  that  can  also  achieve  optimization  if  you  have 
knowledge  of  the  command  sequence  in  the  near  future  [25,27,28]. 
Another  method  used  to  avoid  singularities  is  to  use  null  motion  to 
first  reorient  the  CMGs  to  desired  gimbal  positions  that  are  not  near 
singular  configurations  [29].  By  definition,  null  motion  is  motion  of 
the  CMGs  that  results  in  a  net  zero  torque.  Null  motion  only  exists 
when  more  than  three  SGCMGs  are  used.  The  extra  degrees  of 
freedom  provided  by  a  redundant  configuration  with  more  than  three 
SGCMGs  are  used  to  execute  the  null  motion.  Despite  the  massive 
amount  of  research  done  on  CMGs,  precision  control  with  CMGs 
(especially  only  three  SGCMGs,  a  nonredundant  configuration)  is 
still  an  unsolved  problem  [28,30,31J.  This  research  paper  inves¬ 
tigates  a  modification  of  the  four-CMG  pyramid.  Here,  only  three 
CMGs  are  used;  thus,  the  problem  is  nonredundant:  three-axis 
control  is  accomplished  by  exactly  three  degrees  of  freedom.  There 
will  not  be  an  extra  degree  of  freedom  to  use  a  constraint  equation  for 
singularity  avoidance.  Furthermore,  this  paper  will  describe  methods 
to  exactly  follow  the  commanded  path  to  minimize  tracking  errors 
rather  than  a  path  that  avoids  singularities. 

Typical  CMG  output  torques  are  on  the  order  of  hundreds  to 
thousands  of  times  the  torque  output  of  reaction  wheels,  another  kind 
of  momentum  exchange  attitude  control  actuator.  A  unique  challenge 
of  CMG  implementation  remains  the  mathematical  singularity. 
Arguably,  the  most  common  configuration  for  a  skewed  array  of  four 
CMGs  is  the  pyramid  array  where  the  four  CMGs  are  skewed  at  an 
angle  of  fi  =  54.73  deg,  resulting  in  an  optimal  spherical 
momentum  capability  [j_]  requiring  internal  singularity  avoidance 
[1-8,32-37],  The  desire  is  often  stated  as  an  equivalent  (spherical), 
maximized  momentum  capability  (not  singularity-free)  in  all 

directions  based  on  the  {+  +  ++}  or  { - }  OH  and  AH 

saturation  singularities,  where  all  four  CMGs  are  pointing  in  the 
same  direction  (Fig.  2).  The  typical  design  approach  may  be 


II.  Torque  Generation  and  Singularities 

To  achieve  a  specified  output  torque  from  a  CMG  array,  a 
command  must  be  sent  to  the  gimbal  motor.  Equations  0-4)  derive 
this  relationship  for  i  =  n  CMGs  normalized  by  one  CMG’s  worth  of 
momentum  (1/7).  First,  write  the  angular  momentum  vectors  hx,  hy, 
and  hz  in  x,  y,  and  z  directions  composed  of  components  of  angular 
momentum  contributions  of  the  three  CMGs  combining  to  form  the 
overall  system  angular  momentum  vector  H.  It  is  desirable  to  express 
the  rate  of  change  of  angular  momentum  H  in  order  to  use  the  first- 
principles  Newton-Euler  relationship  between  the  rate  of  change  of 
angular  momentum  and  the  applied  external  torques.  The  rate  of 
change  of  angular  momentum  H  may  be  decomposed  into  the 
gradient  with  respect  to  the  gimbal  angles  9H / 9#,  multiplied  by  the 
rate  of  change  of  the  gimbal  angles  {6 1.  CMGs  are  inclined  so  gimbal 
planes  form  skew  angles  with  respect  to  the  xy  plane  as  depicted  in 
Fig.  2,  where  1=1  —  3  for  three  CMGs..  Begin  by  writing  equations 
for  each  momentum  vector  in  xyz  coordinates  for  three  CMGs 
normalized  by  1 H: 

hx  =  cos  03—  cos  +  cos  /J2  sin  92  t 

hy  =  cos  fi3  sin  d3—  cos  ft  sin  9t—  cos  d2  >  (la) 

hz  =  sin  sin  9t  +  sin  /J2  sin  92  +  sin  fi3  sin  93  ) 


hx  =  cos  93  —  cos  +  cos  /3  sin  02  hy  =  cos  yS(sin  03  —  sin  9{) 

—  cos  92 hz  =  sin  /3(sin  9t  +  sin  92  +  sin  93)  >  (lb) 


H  =  hx.r  +  hyv  +  hz£  (2a) 


{+  +  +} 

HS=3H 


Fig.  2  Individual  CMG  momentum  directions  for  a  3/4  CMG  skewed 
array  in  a  singular  configuration  of  gimbal  angles. 
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dhz 

dd; 

sin  /I  cos  6t 

sin  /J  cos  92 

sin  yS  cos  93  

[A] 

=  [A]  (2b) 


The  [A]  matrix  (containing  gimbal  angles  9t  and  skew  angles  /);) 
must  be  inverted  to  find  the  required  CMG  gimbal  command  for 
commanded  output  torque  per  Eq.  (3).  The  Newton-Euler  relation 
relates  generated  torque  to  the  timed  rate  of  change  of  angular 
momentum  of  the  spacecraft  system.  A  CMG  absorbs  momentum 
change  H,  causing  an  equal  and  opposite  change  in  momentum  on 
the  spacecraft.  For  n  CMGs,  the  general  relation  is 
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9H90 

JeTt 


9H 

lwi 


i  =  l 


[A]-‘{H}  =  [A]-1[A]{0}  =  W 


(3) 


For  some  combinations  of  gimbal  and  skew  angles,  the  [A]  matrix 
columns  can  become  linearly  dependent.  At  these  combinations  of 
skew  and  gimbal  angles,  the  determinant  of  the  [A]  matrix  becomes 
zero,  leading  to  singular  inversion  as  follows: 

det[A]  =  s/3{s92[s(9l  +  92)]  +  c/3c92[s(93  —  93)  +  2c$|C#3c-J6]} 

(4) 


where  s  =  sine,  and  c  =  cosine. 

III.  Three  of  Four  Skewed  Control  Moment 
Gyroscope  Array 

The  3/4  CMG  array  modifies  the  commonly  studied  four-CMG 
skewed  pyramid.  A  minimum  of  three  CMGs  is  required  for  three- 
axis  control,  and  the  fourth  is  often  used  for  singularity  avoidance. 
With  the  3/4  array,  only  three  CMGs  are  used  for  active  attitude 
control  with  the  fourth  CMG  held  in  reserve  for  robust  failure 
properties.  The  fourth  CMG  is  not  active  (dormant),  with  no 
electrical  power  applied  to  its  gimbal  motor  until  required  by  the  loss 
of  a  failing  CMG.  When  a  CMG  fails,  it  may  be  despun  (controlled) 
and  then  depowered  while  the  spare  CMG  is  powered  and  spun  up 
(again  controlled).  Thus,  the  fourth  CMG  substitutes  the  failing 
CMG,  maintaining  a  3/4  CMG  array  configuration.  The  3/4  CMG 
array  remains  the  focus  of  this  research,  and  the  results  here  can  be 
equally  applied  to  the  new  (substituted)  3/4  CMG  array.  Experi¬ 
mental  verification  will  be  provided  in  later  sections  using  a 


spacecraft  testbed  with  a  3/4  CMG  array  containing  a  balance  mass 
in  the  place  of  the  fourth  CMG  (Fig.  2). 

The  approach  taken  by  the  authors  is  to  first  optimize  the  3/4 
skewed  array  (system  of  four  CMGs  in  a  skewed  configuration  where 
only  three  CMGs  operate  for  attitude  control)  geometry  itself  by 
choosing  the  skew  angle  that  provides  the  greatest  singularity-free 
momentum.  At  this  optimal  singularity-free  skew  angle,  the  3/4 
CMG  array  can  operate  at  momentum  values  less  than  the 
singularity-free  threshold  without  any  kind  of  singularity  avoidance 
scheme.  Furthermore,  utilization  of  mixed  skew  angles  can  rotate  the 
work  space  to  maximize  momentum  in  a  preferred  direction,  again 
singularity-free.  Yaw  is  the  preferred  direction  in  this  study.  A  direct 
comparison  with  the  traditional  optimal  spherical  skew  angle  will 
demonstrate  the  dramatic  improvement  in  torque  capability  of  the 
CMG  array.  Analytical  derivation  is  followed  by  heuristic,  geometric 
analysis,  and  then  validation  via  experimentation  on  a  realistic 
spacecraft  simulator  in  ground  tests. 

IY.  Analysis 

Singular  combinations  of  gimbal  angles  and  skew  angles  can  be 
determined  analytically  by  examining  the  determinant  of  the  [A] 
matrix.  Recall 

det[A]  =  s/3{s92[s(9l  +  92)]  +  cfic92[s(83  —  9t)  +  2c93c93cfl]} 

When  the  determinant  goes  to  zero,  the  matrix  has  linearly  dependent 
columns  resulting  in  singular  inversion.  There  are  six  cases  (with 
multiple  subcases)  that  result  in  a  singular  [A]  matrix  (less  than  full 
rank)  with  /),  =  /l,  where  each  singular  case  is  caused  by  a 
component  quantity  being  equal  to  zero  resulting  in  the  total  quantity 
not  being  invertible  (see  Table  j_): 


Case 

1 

2 

3 

4 

5 

6 


Table  1  Six  singular  cases  where  det[A]  =  0 


Value 
sin  P  =  0 

sin(@!  +  93)  =  sin($3  —  03)  =  0 

sin  $2[sin($i  +  03)]  +  cos  /i  cos  (92[sin(03  —  8t )  +  2  cos  9  cos  63  cos  /J]  =  0 
sin  92  =  0 

sin(9[  +  93)  =  cos  92  =  0 
sin(9[  +  83)  =  cos  82  =  0 


z 
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Case  1 : 

sin /I  {sin  02[sin(0!  +  $3)]  +  cos (6 cos $2[sin($3  —  9{) 

0 

+  2  cos  9l  cos  93  cos  P\} 

— >•  sin  P  =  0  (5) 

Case  2: 

sin /3{sin  (^[sin^  +  #3)]  +  cos /J  cos  #2[sin($3  —  9{) 

0  0 

+  2cos  9{  cos  93  cos yS] }  — >•  sini^y  +  93)  =  0  and 
sin(03  —  9\)  =  0  (6) 


before  reaching  a  singular  state  for  three  CMGs.  To  avoid  singu¬ 
larities  and  maintain  attitude  control,  spacecraft  control  torque  is 
limited  to  less  than  0. 15 H.  Thus,  using  the  ubiquitous  skew  angle  of 
54.73  deg,  which  was  optimized  for  spherical  maximum  momen¬ 
tum  (not  singularity  free)  with  four  SGCMGs,  in  the  case  of  three 
CMGs  only,  15%  of  one  CMG’s  torque  is  achievable  singularity  free. 
Using  the  new  skew  angle  of  /l  =  90  deg  (see  Figs.  4  and  5),  which 
was  optimized  for  singularity-free  momentum  space,  1.0 H  (100% 
the  momentum  capability  of  one  CMG)  is  achievable  singularity- 
free  in  any  direction.  Theoretically,  singularity-free  momentum 
is  increased  from  roughly  0.15/7  to  1H:  (IH  —  0.15//)/ 
0.15//  >  500%.  This  theoretical  claim  must  be  validated.  To  do  so. 
Sec.  V  will  provide  simulated  visual  depictions  of  the  singularity 
surfaces  in  the  momentum  space  demonstrating  a  lack  of 
singularities  in  the  momentum  space  less  than  1//,  and  then 
Sec.  VII  will  provide  experiments  that  command  maneuvers 
resulting  in  angular  momentum  trajectories  that  exceed  0. 15 H  while 
monitoring  the  condition  of  the  [A]  matrix  for  potential  singularities. 
If  the  inverse  condition  number  of  the  [A]  matrix  ever  goes  to  zero,  the 


Case  3: 


sin  /l{sin  02[sin(0[  +  03)]  +  cos  p  cos  6b[sin(#3  —  9{)  +  2  cos  9{  cos  93  cos  P\} 

0 

—¥■  sin //2[sin( //j  +  03)]  +  cos/1  cos  $2[sin((93  —  9X)  +  2cos$!  cos  93  cos /l]  =  0  (7) 


Case  4: 

sin/3{  sin  92  [sin(0[  +  #3)]  +  cos  P  cos  $2[sin(#3  —  9{) 

0  0 

+  2  cos  9 1  cos  93  cos  P\  [  — >  sin  92  =  0 
and  cos  P  cos  92  =  0  (8) 

Case  5: 


0  0 


+  2  cos  9{  cos  93  cos  /!]}  — >  sin( +  93)  =  0 
and  cos  P  cos  92  =  0  (9) 

Case  6: 

sin y0{  sin  92  [sin^  +  03)] 

0 

+  cos  P  cos  @2[sin((93  —  93)  +  2  cos  9{  cos  93  cos  P\) 

' - ' 

0 

— >  sin  92  =  0  and  sin(@3  —  9{ ) 

+  2  cos  9X  cos  93  cos  p  =  0  (10) 

There  are  a  few  trivial  cases.  Nontrivial  cases  may  be  analyzed  as 
follows.  In  general,  for  a  given  skew  angle,  each  case  produces 
gimbal  angle  combinations  that  result  in  det[A]  =  0.  These  gimbal 
combinations  may  be  used  to  calculate  the  resultant  momentum  at  the 
singular  condition.  Minimum  singular  momentum  values  may  then 
be  plotted  for  iterated  skew  angles  0  deg  <P  <  90  deg.  Having 
established  the  minimum  value  of  momentum  at  singular  combi¬ 
nations  of  gimbal  angles,  it  is  not  possible  to  become  singular  at 
momentum  less  than  these  values  (Fig.  3).  Thus,  the  result  is  the 
maximum  singularity-free  momentum  space.  As  we  see  in  Fig.  3,  the 
extensively  studied  skew  angle  (P  =  54.73  deg),  which  is  originally 
optimized  for  maximal  momentum  performance  of  four  CMGs  (not 
singularity-free  momentum),  yields  roughly  0.15//  momentum 


CMGs  have  encountered  a  singularity.  It  will  be  shown  that  the 
CMGs  remain  singularity  free  as  claimed  (see  Fig.  4). 


V.  Heuristics 

The  preceding  analysis  reveals  singularity-free  operations  less 
than  1H  in  all  directions  by  implication.  While  useful,  the  analysis 
certainly  does  not  yield  much  intuition  for  the  attitude  control 
engineer  to  design  safe  momentum  trajectories  through  the  momen¬ 
tum  space.  Are  there  directions  that  can  exceed  1 H  singularity-free? 

Advances  in  computer  processor  speeds  make  a  heuristic  approach 
readily  available.  Consider  rotating  a  vector  360  deg  creating  a  CMG 
gimbal  cutting  plane  (discretized  at  some  interval).  Then,  rotate  the 
gimbal  plane  360  deg,  creating  a  lattice  of  discrete  points  forming  a 
solid,  filled  sphere.  This  lattice  provides  discretized  points  to  analyze 
CMG  array  momentum.  This  may  be  done  in  embedded  loops  of 
computer  code.  Each  discrete  point  corresponds  to  a  set  of  three 


Fig.  4  Skew  angle  optimization  for  maximum  singularity-free 
normalized  momentum. 
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Fig.  5  3/4  CMG  skewed  array  momentum  cutting  planes  with 
/}  =  90  deg. 


Fig.  6  Singular  surfaces  for  a  traditional  3/4  CMG  skewed  array  at  a 
skew  angle  of  54.73  deg.. 


Fig.  7  Skew  iteration  at  70  deg  (top),  80  deg  (lower  left),  and  90  deg 
(lower  right). 


Fig.  8  Mixed  skew  angles  (see  Table  2)  highlighting  case  2 
{90,0,90  degjC 


coordinates  or,  equivalently,  three  gimbal  angles.  At  each  discrete 
point,  the  singularity/nonsingularity  of  [A]  is  established.  At  singular 
points,  the  normalized  magnitude  of  angular  momentum  may  be 
calculated.  A  point  may  be  plotted  at  the  magnitude  of  the  momen¬ 
tum  in  the  singular  direction.  This  results  in  a  three-dimensional 
singularity  map  granting  intuition  for  maneuvering  in  the  momentum 
space  (e.g..  Figs.  6-8).  When  the  singular  condition  exists,  one,  two, 
or  three  CMGs’  angular  momentums  may  be  generating  momentum 
in  the  same  direction,  and  this  fact  leads  to  the  nomenclature  listed  in 
Fig.  5,  which  is  ubiquitously  used  in  the  literature  (e.g.,  [37]).  If  all 
three  CMGs  point  in  the  same  positive  or  negative  direction  (OH  and 
3 H),  this  represents  the  maximum  capability  of  the  array,  and  the 
singularity  is  referred  to  as  saturation.  Other  singular  cases  that  occur 
when  one  CMG  is  pointing  opposed  to  the  other  two  CMGs  (1/7)  are 
referred  to  as  internal. 

Singular  surfaces  plotted  in  the  momentum  space  (e.g..  Figs.  6-8) 
result  from  several  kinds  of  singularities:  OH,  1 H,  2 H,  and  3 H, 
depicted  in  Fig.  5.  Typically,  the  optimal  spherical  skew  angle  is 
determined  by  the  outer,  saturation  singular  surface  alone  (OH  and 
3 H  in  the  case  of  a  3/4  CMG  array).  To  find  the  maximum 
singularity-free  skew  angle,  the  analysis  must  also  account  for  the  1 H 
and  2 H  internal  singularities. 

It  may  be  noted  here  that  the  often  used  optimal  spherical  skew 
angle  (fi  =  54.73  deg)  that  results  in  equivalent  momentum  in  all 

directions  is  derived  using  the  OH  { - }  and  4 H  {+  +  ++} 

singular  surfaces  of  a  four-CMG  skewed  array.  The  OH  and  AH 
singular  surfaces  are  the  saturation  surfaces  that  result  from  all  four 
CMGs  pointing  in  the  same  direction,  and  their  corresponding 
saturation  surfaces  for  the  3/4  array  are  the  OH  and  3 H  singular 
surfaces.  In  the  case  of  a  3/4  CMG  array,  the  singular  surfaces  are 
depicted  in  Fig.  6  for  the  often  used  optimal  spherical  momentum 
skew  angle  (yS  =  54.73  deg). 

Note  that  the  spherical  nature  of  the  external,  saturation  singular 

surfaces  {+  +  +}and{ - }  typified  by  the  ubiquitous  four-CMG 

pyramid  array  skewed  at  /)  =  54.73  deg  (Fig.  1)  are  nearly 
maintained  in  the  3/4  array.  The  internal  1 H  and  2 H  singularities  are 
quite  a  problem,  since  they  occupy  a  large  portion  of  the  momentum 
space.  The  combined  singularity  hypersurface  makes  it  difficult  to 
see  a  clear  momentum  path  around  the  origin. 

The  typical  design  methodology  might  use  the  shape  of  the  3 H 

outer  momentum  surfaces  {+  +  +}  and  { - }  to  define  /3  = 

54.73  deg  as  the  optimal  spherical  momentum  skew  angle.  The 
attitude  control  engineer  would  be  left  with  the  task  of  maneuvering 
in  this  crowded  momentum  space  while  trying  to  avoid  any  point  on 
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the  singular  hypersurface.  Striking  a  singular  point  results  in  (at  least 
temporary)  loss  of  attitude  control.  Also  note  the  maximum 
momentum  capability  is  less  than  3 H.  Per  Fig.  3,  the  maximum 
singularity-free  momentum  capability  using  fi  =  54.73  deg  is 
0.154868//.  If  the  skew  angle  were  increased  to  90  deg.  Fig.  3  skew 
angle  optimization  indicates  that,  theoretically,  singularity-free 
momentum  would  be  increased  1.0 H.  Consider  the  singular  hyper¬ 
surfaces  for  heuristic,  geometric  observations. 

It  was  established  in  Sec.  IV  that  a  skew  angle  of  90  deg 
theoretically  results  in  a  singularity-free  momentum  space  of  IH  for 
the  3/4  CMG  array.  Repeating  the  numerical,  heuristic  singularity 
analysis  as  the  skew  angle  increases  is  very  revealing.  Figure  7 
displays  the  singular  momentum  space  for  a  skew  angle  of  70,  then 
80,  then  90  deg.  Similar  plots  incremented  every  fi  =  5  deg  are 
contained  in  [38],  while  these  three  examples  in  Fig.  7  are  sufficient 
to  demonstrate  the  behavior  of  increasing  the  skew  angle.  Notice  how 
the  1 H  and  2 H  singular  surfaces  move  away  from  the  origin  and 
gradually  converge  into  each  other  as  the  skew  angle  increases  to 
90  deg.  Fewer  singular  surfaces  are  obviously  beneficial,  but  the 
vacancy  of  the  center  of  the  momentum  space  is  especially 
significant,  since  it  creates  a  significant  portion  of  the  momentum 
space  that  is  free  of  singular  conditions  (free  for  maneuver  without 
regard  to  singularities).  Furthermore,  it  will  be  shown  in  the  next 
section  that  the  large  area  of  singularity-free  momentum  space  can  be 
rotated  via  mixed  skew  angles  to  emphasize  a  preferred  axis  of 
rotation.  The  yaw  axis  is  the  preferred  axis  of  maneuver  in  this  study, 
but  notice  in  Fig.  7  that  greater  angular  momentum  (2 H)  capability 
exists  in  the  roll  direction  for  this  geometry  of  CMG  array,  while 
pitch  maneuvers  can  achieve  1 H  singularity  free.  It  is  desirable  to 
rotate  the  direction  of  maximum  singularity-free  momentum  space  to 
the  yaw  direction.  It  will  be  seen  that  mixed  skew  angles  can  rotate 
the  momentum  space  to  emphasize  roll  or  pitch. 


VI.  Momentum  Space  Rotation:  Mixed  Skew  Angles 

Typically,  skewed  CMG  arrays  use  identical  skew  angles  for  each 
CMG  (Pj  =  ft).  By  using  mixed  skew  angles,  the  singularity-free 
football- shaped  space  can  be  reoriented  to  place  the  2 H  momentum 
capability  in  the  yaw  direction.  Six  possible  momentum 
reorientations  are  possible  by  laying  down  momentum  planes  from 
90  to  0  deg,  as  listed  in  Table  2,  resulting  in  rotations  of  the  momen¬ 
tum  space  depicted,  respectively,  in  Fig.  8.  This  section  refers  to  fixed 
geometries  and  does  not  investigate  the  adaptive  skew  angle  [39], 
especially  since  the  adaptive  skew  angle  provides  an  extra  degree  of 
freedom,  while  this  research  seeks  to  investigate  what  can  be  done 
with  nonredundant  configurations. 

Notice  (in  Fig.  8)  that  three  options  for  mixed  skew  angles  result  in 
the  original  momentum  space  rotated  about  z  such  that  x  -o-  y,  while 
two  other  options  generate  spherical  momentum  space  filled  with 
significant  internal  singularities.  Notice  the  center  of  the  momentum 
space  is  clogged  with  singular  surfaces  such  that  those  two  figures  are 
blackened  in  the  center.  Our  difficulty  seeing  the  center  is  indicative 
of  difficulties  steering  a  momentum  vector  through  that  space 
without  striking  a  singular  surface  for  those  two  mixed  skew  angle 


Table  2  Six  possible  combinations  of  mixed 
skew  angles  laying  one  or  two  momentum 
cutting  planes  from  0  to  90  deg- 


Pi’  deg 

Pi,  deg 

Pi’ deg 

0 

90 

90 

90S 

os 

90S 

90 

90 

0 

0 

0 

90 

0 

90 

0 

90 

0 

0 

“Corresponding  singular  hypersurfaces  are  depicted  in 
Fig.  8,  starting  from  the  upper  left  and  going  to  the 
lower  right 

bSpecial  attention  is  brought  to  the  case  of 
{ft,. fi2, ft,}  T  =  {90,0,90  deg}T. 


Z 


Fig.  9  {/?!,  /t2,  Pi)  =  {90, 0, 90  deg)  mixed  skew  angle  momentum 

cutting  planes  depicted  for  a  3/4  CMG  skewed  array. 


combinations.  This  is  an  undesirable  feature  similar  to  the  original 
case  of  four  CMGs  skewed  at  54.73  deg  (Fig.  _1)  or  three  CMGs 
skewed  at  54.73  deg  (Fig.  6)  where  only  0.15//  was  available 
singularity-free.  This  study  seeks  geometries  with  large  areas  of  the 
momentum  space  that  are  free  of  singularities. 

One  successful  reorientation  /32,  P3}  =  {90,  0, 90  deg}  is 
accomplished  by  simply  sliding  the  second  CMG  from  90  to  0  deg, 
resulting  in  singularity-free  yaw  momentum  increase  with  a 
favorably  vacant  internal  singularity  structure.  Figure  9  depicts  a  3/4 
CMG  array  skewed  at  mixed  skew  angles  yS  (  =  90  deg,  /t2  =  0  deg, 
and  p3  =  90  deg.  Trajectories  originating  at  the  origin  of  the 
momentum  space  have  1 H  spherical  momentum  capability  and  2 H 
momentum  capability  about  yaw  (z)  singularity-free.  Theoretically, 
singularity-free  momentum  is  increased  to  2 H  from  0.15//  compared 
with  the  case  of  three  CMGs  skewed  at  54.73  deg  (Figs.  3  and  7): 
(2H  -0.15//)/ 

0.15//  >  1000%.  This  claim  must  be  experimentally  validated. 
Section  VIII  provides  experiments  that  command  maneuvers 
resulting  in  angular  momentum  trajectories  that  exceed  0.15//  and 
even  exceed  1 H  (the  demonstrated  limit  of  singularity-free  momen¬ 
tum  space  with  nonmixed  skew  angles).  Note  that  momentum 
trajectories  that  are  initiated  from  points  near  (0,0,2)  can  traverse  to 
(0,0, —2),  resulting  in  —4 H  being  stored  in  the  CMG  array 
producing  +4 H  momentum  imparted  to  the  spacecraft  about  yaw 
singularity  free. 

VII.  Experimental  Verification: 

Skew  Angle  Optimization 

Experimental  verification  of  the  singularity-free  skew  angle  is 
performed  on  a  free-floating  three-axis  spacecraft  simulator  to 
demonstrate  singularity-free  operations.  The  free-floating  spacecraft 
testbed  is  referred  to  as  TASS2  (three-axis  spacecraft  simulator) 
indicating  its  heritage  as  the  second  such  testbed  developed  at  the 
Naval  Postgraduate  School.  The  first  testbed,  TASS  1  [40],  was  much 
smaller  and  did  not  use  CMGs.  A  description  of  the  experimental 
hardware  is  based  on  [38],  for  which  the  graphics  were  directly 
replicated  here  (Fig.  JB).  Figure  JB  shows  the  actual  picture  of  the 
three-axis  simulator  named  as  TASS2.  The  spacecraft  bus  simulator 
is  supported  on  a  spherical  air  bearing  to  simulate  a  weightless 
environment.  A  thin  film  of  compressed  air  is  injected  between  a 
spherical  ball  and  a  mating  spherical  cup.  This  thin  film  of  air  creates 
an  essentially  frictionless  lubrication  layer  between  the  ball  and  cup. 
When  test  articles  mounted  to  the  ball  segment  are  balanced,  such 
that  their  aggregate  center  of  gravity  corresponds  with  the  center  of 
rotation  of  the  ball,  rotational  motions  of  the  ball  and  test  articles 
match  those  of  an  object  with  similar  inertial  properties  free  falling 
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Fig.  10  TASS2,  where  IMU  refers  to  inertial  measurement  unit. 


though  space.  Because  it  is  important  to  minimize  the  disturbance 
torque  from  any  imbalance,  three  servolinear  stages  with  lumped 
masses  may  be  employed  for  automatic  mass  balancing.  The  three 
linear  stages  are  placed  parallel  to  the  three  spacecraft  body  axes  but 
are  not  used  in  this  study. 

The  spacecraft  testbed  payload  is  an  optical  relay  for  laser  beams. 
The  lower  deck  of  the  spacecraft  is  a  bench  model  spacecraft  with 
CMGs,  active  mass  balancing  system,  fiber  optic  gyroscopes  for 
attitude  rate  sensing,  and  typical  onboard  computers.  The  upper  deck 
of  the  testbed  is  an  experimental  laser  optical  relay  designed  to  accept 
a  laser  from  the  ground,  aircraft,  or  space-based  source  and  relay  the 
laser  to  an  uncooperative  target  on  the  ground,  in  the  air,  or  in  space. 
Description  of  the  mission  and  some  details  about  experimental 
hardware  may  be  found  in  [40,41J.  This  laser  relay  mission  demands 
rapid  target  acquisition  yaw  maneuvers  in  minimum  time  followed 
by  fine  pointing  during  target  tracking. 

To  first  demonstrate  singularity-free  operations  (verify  the  large 
maneuver  space  is  free  of  singular  conditions  as  theoretically 
predicted  earlier  in  this  paper),  experimental  maneuvers  were  per¬ 
formed  with  a  +5  deg  yaw  maneuver  in  4  s  followed  by  a  —5  deg 
yaw  maneuver  in  4  s.  The  attitude  is  then  regulated  to  zero  while  the 
CMG  continues  to  output  torque  to  counter  gravity  gradient 
disturbances  typical  of  imbalanced  ground  test  spacecraft  simulators. 
The  testbed  has  an  autobalancing  device  to  eliminate  this  gravity 
disturbance  torque,  but  the  device  was  disengaged  to  permit  the 
experiment  to  explore  more  of  the  momentum  space  and  to  simplify 
assertions  of  CMG  performance,  since  CMGs  will  be  the  only  torque 
actuators.  During  previous  research,  this  maneuver  (+5  deg  yaw 
maneuver  in  4  s  is  followed  by  a  —5  deg  yaw  maneuver  in  4  s)  was 
performed  using  a  CMG  geometry  with  skew  angle  of  57  deg  (not 
depicted).  The  CMG  array  became  singular,  and  the  testbed 
spacecraft  attitude  control  was  lost,  motivating  this  current  study. 
The  skew  angle  was  increased  to  90  deg  for  all  three  CMGs,  and  the 


Fig.  11  Experiment:  +5  deg  yaw  in  4  s,  —5  deg  yaw  in  4  s,  and  then 
regulate  for  5  min. 


time  (sec) 

Fig.  12  Experiment:  +5  deg  yaw  in  4  s,  —5  deg  yaw  in  4  s,  and  then 
regulate  for  5  min.  Skew  angles  equal  90  deg. 


identical  experiment  was  repeated.  Notice  in  Figs.  IT  and  12  that  the 
maneuver  is  performed  and  the  testbed  is  regulated  for  5  min.  without 
striking  any  singular  surfaces.  Momentum  magnitude  and  the  inverse 
of  the  condition  of  the  [A]  matrix  verify  this  assertion. 

Notice  what  happens  when  the  same  momentum  trajectory  is 
placed  in  the  context  of  the  theoretical  singular  momentum  space  of 
the  optimal  spherical  skew  angle  /)  =  54.73  deg  (Fig.  _13).  The 
internal  singular  surfaces  are  depicted  individually  for  ease  of 
visualization.  This  momentum  trajectory  is  constantly  close  to  the 
internal  singular  surfaces  and  quickly  strikes  a  singular  surface.  A 
corresponding  singular  surface  exists  for  a  skew  angle  of  57  deg. 
Prior  experiments  using  f)  =  51  deg  went  singular  and  resulted  in 
the  loss  of  attitude  control  when  the  momentum  trajectory  struck  this 
corresponding  singular  surface. 


VIII.  Experimental  Verification: 

Momentum  Space  Rotation 

Next,  experiments  were  performed  with  mixed  skew  angles  to 
orient  the  maximum  singularity-free  momentum  capability  about  the 
yaw  axis,  as  seen  in  Figs.  ]_4  and  15.  The  first  experiment  (Fig.  1_1) 
verified  that  momentum  trajectories  could  exceed  1H  in  roll,  but  care 
was  taken  not  to  exceed  1 H  in  yaw,  since  1 H  defines  the  yaw  limit  for 
singularity-free  operations  of  the  {90, 90, 90  deg}  configuration.  To 
verify  momentum  space  rotation  here,  we  seek  to  verify  that  we  can 
exceed  1 H  about  yaw,  so  the  yaw  maneuvers  were  increased  1 60%  in 
the  same  duration  from  5  to  13  deg  in  only  4  s.  This  demands 


Fig.  13  Fictional  experiment:  Fig.  IT  placed  in  context  of  54.73  deg 
momentum  envelope. 
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Fig.  14  Experimental  results:  +13  deg  yaw  in  4  s,  and  —13  deg  yaw 
in  4  s  performed  in  {/?i,  j82, /I3}  =  {90,0,90  deg}  mixed  skew  angle 
configuration,  the  momentum  trajectory  is  placed  in  the  context  of 
theoretical  singular  hypersurface. 


Fig.  15  Experimental  results:  +13  deg  yaw  in  4  s,  and  —13  deg  yaw 
in  4  s  performed  in  {/?!, /?2, /?3)  =  {90,0,90  deg}  mixed  skew  angle 
configuration  (with  maneuver  momentum  and  inverse  condition  of  [A]). 

significantly  more  momentum  change,  specifically  about  yaw.  Per 
Figs.  14  and  15,  the  momentum  is  achieved  singularity-free  and  the 
maneuver  is  performed  without  incident;  thus,  we  see  that  the 
momentum  space  can  be  rotated  to  increase  yaw  maneuver 
capability. 


IX.  Analysis  of  Singularity  Reduction:  Decoupled 
Steering  Commands  Technique 

In  this  section,  we  derive  a  strategy  dubbed  “decoupled  control 
steering”,  where  we  take  advantage  of  the  simplifications  that  arise 
from  the  optimum  singularity-free  skew  angle  /J  =  90  deg,  seeking 
to  yield  maximal  (3 H)  momentum  capability  about  yaw  without  loss 
of  full  attitude  control.  Substituting  the  [A]  matrix  with  f)  =  90  deg 
into  Eq.  (2)  yields 
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Note  that  the  y-momentum  change  equation  has  become 
decoupled  from  the  jc  and  z  equations.  Pitch  momentum  is 


determined  completely  by  gimbal  2.  The  pitch  equation  may  be 
separated  from  the  matrix  system  of  equations.  The  benefit  is  the 
elimination  of  singular  gimbal  commands  for  CMGs  that  are  not  in 
geometrically  singular  gimbal  angle  positions.  Consider  what 
happens  if  the  first  and  third  CMGs  enter  a  combination  of  gimbal 
angles  that  satisfy  cos  8{  sin  93  +  sin  0l  cos  03  =  0.  This  would  not 
result  in  singular  commands  to  CMG  gimbal  2.  CMG  gimbal  2  would 
receive  the  following  command:  02  =  ( 1/  sin  02)hy.  The  individual 
equations  for  each  of  the  three  CMGs  to  be  implemented  are 
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During  singular  conditions,  nonsingular  CMGs  would  operate 
normally  per  their  decoupled  steering  logic.  Take  special  care  not  to 
implement  the  equations  listed  below  in  a  seemingly  equivalent 
matrix  form  [Eq.  (12)].  By  calculating  each  CMG’s  proper  command 
individually  (rather  than  in  a  matrix),  decoupled  control  steering 
allows  nonsingular  CMGs  to  be  properly  commanded  during  periods 
that  other  CMGs  are  singular.  This  can  prevent  loss  of  full  attitude 
control  and  permits  the  momentum  trajectory  to  pass  the  internal 
singular  state. 


X.  Simulation  of  Singularity  Reduction 

Yaw  maneuvers  were  simulated  using  typical  coupled  control  and 
compared  with  the  proposed  decoupled  control  strategy.  First,  a 
+50  deg  yaw  maneuver  is  followed  immediately  by  a— 50  deg  yaw 
maneuver,  and  then  regulation  at  zero.  The  results  of  both  methods 
are  displayed  in  Fig.  1_6.  Notice  the  coupled  implementation  of  the 
Moore-Penrose  pseudoinverse  results  in  large  unintended  roll  each 
time  the  momentum  trajectory  strikes  the  singular  surface. 

On  the  contrary,  notice  how  decoupled  control  smoothly  traverses 
the  singularity  surface  with  negligible  roll  or  pitch  errors.  Since 
analysis  and  simulation  both  indicate  the  proposed  decoupled  control 
technique  should  work,  experimental  verification  was  performed  on 
a  free-floating  spacecraft  simulator  (Figs.  17-19). 

Notice  in  the  simulation  depicted  in  Fig.  JA  that  identical 
simulated  experiments  are  compared,  except  one  of  the  experiments 
used  a  coupled-matrix  steering  law  (dotted  line),  while  the  second 
used  uncoupled,  individual  steering  equations  to  command  each 
CMG  (solid  line).  The  thin  line  is  hard  to  see  when  both  experiments 
travel  identical  paths  through  the  momentum  space,  but  when  the 
coupled-matrix  steering  law  encounters  a  singularity,  loss  of  attitude 
control  results  in  large,  undesired  roll  maneuvers.  Underneath,  it  is 
easier  to  see  the  think  line  that  represents  the  simulation  using  the 
uncoupled,  individual  steering  equations,  which  cleanly  passes 
through  the  singularity  without  loss  of  attitude  control. 


XI.  Experimental  Verification:  Singularity  Reduction 
via  Decoupled  Steering  Technique 

Experiments  were  performed  with  decoupled  control  to  maximum 
momentum  capability  about  the  yaw  axis.  First  note  that  Fig.  T7 
displays  the  ability  of  decoupled  control  steering  to  penetrate  the 
singular  surface  associated  with  the  coupled  [A]  matrix  of  the  CMG 
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Fig.  16  Simulation:  comparison  of  typical  coupled  control  (dotted  line) 
with  decoupled  control  (solid  line).  Notice  roll  ( x )  errors  when  the 
momentum  trajectory  passes  through  the  singularity  resulting  from  loss 
of  attitude  control. 


Fig.  17  Experimentation:  demonstrated  ability  to  pass  cleanly  through 
the  singular  surface  at  1 H  using  proposed  decoupled  control. 

gimbal  angles  and  skew  angle.  This  attribute  is  exploited  with  an 
aggressive  yaw  maneuver  (Fig.  _18).  The  commanded  yaw  maneuver 
angle  was  increased  700%  from  ±5  deg  in  4  s  to  ±35  deg  in  10  s. 
This  demands  significantly  more  momentum  change,  specifically 


0  50  0  10  20 

Time  (s)  Time  (s) 

Fig.  18  Experiment:  +35  deg  yaw  in  10  s,  and  —35  deg  yaw  in  10  s 
with  decoupled  steering. 


Fig.  19  Experiment:  +35  deg  yaw  in  10  s,  and  —35  deg  yaw  in  10  s 
with  decoupled  steering.  CMG  gimbal  angles  and  1  / cond|.t  |. 


about  yaw.  Figure  J_8  displays  the  required  maneuver  is  achieved 
without  incident.  Notice  that  the  coupled  [A]  matrix  was  singular 
twice  during  this  maneuver. 

Typical  coupled  control  steering  would  have  resulted  in  loss  of 
spacecraft  attitude  control.  Instead,  with  decoupled  steering,  you  will 
notice  a  nice  maneuver  despite  the  singular  [A]  matrix.  Attitude 
control  is  not  lost  at  any  time.  Also  notice  the  extremely  high 
magnitude  of  momentum  achieved  without  loss  of  attitude  control 
associated  with  the  passage  of  singularities. 

XII.  Conclusions 

These  experiments  validate  the  much  desired  goal  of  CMG 
attitude  control:  extremely  high  torque  without  loss  of  attitude 
control  associated  with  mathematical  singularity.  The  optimized 
geometry  is  shown  to  increase  singularity-free  torque  capability 
significantly.  The  maximized  singularity-free  momentum  may  be 
rotated  to  a  preferred  predominant  maneuver  direction  by  using 
mixed  skew  angles.  Using  a  proposed  decoupled  control  strategy, 
further  singularity  reduction  is  achieved  that  is  shown  to  allow 
momentum  trajectories  to  cleanly  pass  through  singular  surfaces 
without  loss  of  attitude  control.  These  claims  were  introduced 
analytically,  and  promising  simulations  were  provided.  Finally, 
experimental  verification  was  performed,  demonstrating  maximal 
yaw  maneuvers  that  passed  through  singular  surfaces  that  would 
render  loss  of  attitude  control  using  typical  coupled  control 
techniques. 
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